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Abstract. Resonant particle production, along with many other physical processes 
which change the effective equation of state (EOS) during inflation, introduces a feature 
in the primordial power spectrum which in many models has a step-like shape. We 
calculate observational constraints on resonant particle production, parametrized in 
the form of an effective step height, N e ff and location in fc-space, fcbrcak- Combining 
data from the cosmic microwave background and the 2dF Galaxy Redshift Survey 
yields strong constraints in some regions of parameter space, although the range in 
fc-space which can be probed is restricted to k ~ 0.001 — O.l/i Mpc -1 . We also discuss 
the implications of our findings for general models which change the effective EOS 
during inflation. 
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1. Introduction 

The recent measurements of the temperature anisotropics in the cosmic microwave 
background (CMB) from the Wilkinson Microwave Anisotropy Probe (WMAP) pQ 
[Hj have set a new standard for precision in cosmology. In addition to pinning down 
several key cosmological parameters with unprecedented accuracy, the WMAP data 
have confirmed important aspects of the inflationary paradigm. Quantum fluctuations 
during the inflationary phase are commonly thought to be the seeds of the anisotropics 
in the CMB and large-scale structure in the Universe. The observations are consistent 
with the inflationary mechanism for generating curvature fluctuations on superhorizon 
scales, and with the fluctuations being adiabatic and Gaussian as in the simplest single- 
field models ||. The fluctuations are conveniently characterized by the power spectrum 
of the curvature perturbations Vn(k), and in the simplest models it is proportional to 
some power of the comoving wavenumber k of a given Fourier mode, Vn{k) oc k ns " 1 , 
where n s is the so-called scalar spectral index. For n s = 1 one has the case of scale- 
invariant fluctuations. The data favour a slight red tilt of the scalar spectral index, 
n s < 1, again consistent single-field models [2j. 

However, there are several possible extensions to the simplest models for inflation, 
and in most of them deviations from the simple power-law form of the primordial 
power spectrum are produced. Models which produce features in the primordial power 
spectrum include variants of extended inflation jjj, models with multiple episodes of 
inflation jHj, and phase transitions during the inflationary phase [H]. In this paper we 
consider another proposal, resonant particle production during inflation, which extends 
the single-field models by coupling the inflaton to a massive field (see e.g. [TOD- By this 
mechanism energy can be extracted from the inflaton field during inflation, altering its 
classical motion and producing features in the primordial power spectrum. 

Our aim in this paper is to calculate in detail the shape of the primordial power 
spectrum in the presence of resonant particle production, and to constrain this class of 
models using a collection of the latest cosmological data, including the WMAP CMB 
power spectrum and the 2dF Galaxy Redshift Survey (2dFGRS) galaxy power spectrum 
[TT] . It is important to bear in mind that V-ji(k) is not directly observable. Whether we 
look at the CMB or the large-scale distribution of matter, what we see is the primordial 
fluctuations modified by astrophysical processes. Thus, our ability to constrain the 
primordial fluctuations depends on our understanding of the physics involved in the 
processing, and on our knowledge of the relevant parameters controlling the physics, for 
instance the matter density. Combining different cosmological observations is therefore 
essential, since different observations probe different combinations of cosmological 
parameters, and thus combining them breaks parameter degeneracies and allows tighter 
constraints to be obtained. 

The structure of this paper is as follows. In section 2 we describe the model 
and the calculation of the primordial power spectrum in the presence of resonant 
particle production, and in section 3 we comment on the generality of our results. The 
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observations we use, CMB data and the galaxy power spectrum, are described in section 
4, and section 5 discusses the likelihood analysis and the results, before we conclude in 
section 6. 

2. Resonant particle production 

We will work with a simple single field model for the inflaton. We do this because, as 
will be seen, the main interest here is not the slow change in the power spectrum during 
the course of inflation, but rather, what happens if a resonant phenomenon — namely a 
sudden burst of particle production — radically changes the dynamics for a short period 
of time. 

Our specific toy model is that of a single scalar field <ft, dubbed the inflaton, which is 
rolling slowly down its potential in accordance with the slow-roll conditions, dominating 
the energy-momentum tensor. The inflaton is coupled through the mass to a massive 
boson p in such a way that at a certain moment to f becomes massless. Initially p, being 
heavy compared to the natural mass scale set by the expansion rate rn 2 naturoX ~ H 2 , is 
in its vacuum state. The potential is 

V(<f>, (p) = V{4>) + \M{m - (70) V 2 , (1) 

where we will allow for the case where ip represents some degenerate species and we 
denote the number of different species by M as done in [TU] . The phenomenon is strongly 
resonant essentially because any particles created will be red shifted away exponentially 
fast. 

To find the backreaction we will treat the inflaton at the classical level, while the 
p field is fully quantized. Computing the backreaction implies finding the expectation 
value of the energy-momentum tensor associated with the p field. This is, due to 
problems with proper regularization and renormalization, a highly nontrivial problem 
(see e.g. [12] [E] and [H3D-6J for the fermionic case), but as a lowest approximation we 
can use the Hartree approximation. 

2.1. Formalism and equations of motion 

We assume a spatially flat Friedmann-Robert son- Walker (FRW) model overlayed with 
linear scalar perturbations and the metric takes the form [U| 

ds 2 = -(1 + 2A)dt 2 + 2a(t)d i Bdx i dt + a 2 (t)[(l - 2^)5^ + 2d i d j E}dx i dx j , (2) 

where a(t) is the scale factor and A, B, ip and E describe the possible scalar metric 
perturbations. Likewise given iV scalar fields x 1 1 I = 1? • • • > N, they can be split into a 
homogeneous part and a fluctuation 

X I {t,x) = X I {t) + 5 X I {t,x). (3) 

Because of the inherent gauge freedom related to the arbitrariness of coordinates, 
we can eliminate two out of the four metric perturbations by choosing the spatial 



Observational constraints on particle production during inflation 



4 



and temporal hypersurfaces in a proper way [TBI HH1- The physical content is of 
course not changed by such a coordinate — or gauge — transformation. Certain algebraic 
combinations of the metric and matter perturbations in any gauge can be found that 
correspond to a certain perturbation in a specific gauge. These combinations are dubbed 
gauge invariant variables. 

In this paper we will use the spatially flat gauge to describe the field perturbations. 
These are called the Sasaki-Mukhanov variables: 

Q 1 = Sx 1 + ^V- (4) 



Sx 1 + ty. 

The fourier transform Q l k of Q 1 has the evolution equation [2U] 



(5) 



where M. is the effective mass matrix 



M 1 , 



d 2 V 
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(6) 
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A popular measure of the curvature perturbation, which later directly relates to the 
temperature perturbation in the CMBR, is 1Z — the spatial curvature seen by observers 
who are comoving with the energy content. It is conveniently related to the Sasaki- 
Mukhanov variables as [2JJ 

(7) 



n 



i 



X 1 



Xjx j x j 



Qi 



In our case {x 1 } — {0> V 9 }- 

The evolution of the background parameters, taking into account the backreaction 
of the ip field, but not rescattering of the inflation field itself is found to be: 

dV{<f>) 



H 



H 



-3H(j) 

Atx 
3m 2 pl 

47T 



d<j) 



+ gM{m~g(t))^ 2 ), 



rn 



pi 



<f) 2 + 2V(<f>) + AT(v 2 ) + M(m - g<j)) 2 (v 2 ) + ^(V(p 2 ) 



(0 2 ) + ^(v^ 2 ) 



(8) 
(9) 
(10) 



The expectation values (•) will be computed below. 



2.2. Quantization of the (p field 

When we quantize the (p field we will disregard any effects arising from the metric 
perturbations. For example, we will consider 5ip instead of Q\ effectively evolving 
according to equation © without the second term in equation @. This is reasonable 
as long as we are looking at wavelengths shorter than the horizon size. We will see 
below that indeed the bulk of the energy is deposited in sub-horizon modes justifying 
the approach. The backreaction of 0-modes on the yj-field can also be disregarded, 
because it is an almost homogeneous light field. 
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When quantizing ip we have chosen to use the rescaled field a 3 ^ 2 ip. This is done 
in order to get rid of trivial prefactors arising from the expansion of background 
spacetime. Working in the Heisenberg picture the (p field can be canonically quantized 
and decomposed in terms of annihilation and creation operators 

a 3/2 = — ^ J d 3 M fc (a 3 /V) (t)e-^ + a£ (a 3 / 2 ^)* (t)e +ik *, (11) 

satisfying the usual commutation relations 

[a k , a k >] = [%, a\,\ = 0, [a k , %,} = 5 3 (k - k'). (12) 

By construction (p k has to be a solution to the classical Klein-Gordon equation of a 
massive field with time-dependent mass 

(a 3/ V) " + u 2 k (a 3 / V) = 0, (13) 

u? h= # + ( m - g <j>?-\lP-lH. (14) 

Except for the very small time interval in which the particles are produced, the field is 
heavy compared to the 'Hubble mass' 

(m - g4>) 2 » \H\ (15) 

and the vacuum fluctuations are exponentially damped on large scales. When (m — 
g<p) 2 < \H 2 particle production occurs, and we can assume it to be the dominant 
effect. Therefore we will disregard vacuum polarization altogether. To understand the 
qualitative behaviour of the dynamics, we do not need to include the terms coming from 
the metric equation © nor terms in equation (JT4*j) that are dependent on the Hubble 
parameter. 

Initially (p k is in the vacuum state of a heavy bosonic field 



a" v/v = ^i=exp (-i J u} k dt^J , (16) 



7 3/2 _ 

while after some moment to, where 



m - g<f) = 0, (17) 

particle production, because of interaction with the inflaton, has occurred. Generally 
ip k can be written in terms of 'adiabatic' solutions to equation (JT3j) 



a 3/2( Pk(t) = "^§= ex P ( -i / ^kdt^j + J^— exp ^+i J LU k dt^j . (18) 

Then initially a k — 1, /3 k — 0, while at later times the effective mass (m — g<p) only 
changes slowly and equation (fTHj) is a solution to equation (fT3*j) with constant a k , (3 k . 
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2. 3. Characteristics of the (p field 

Choosing a special boundary condition for equation (fT3|) is equivalent to singling out 
a special set of observers (and annihilation and creation operators). We have already 
mentioned one, namely equation (j!6|) . which is related to the vacuum seen by observers 
in the infinite past. Another interesting vacuum is the co moving vacuum related to the 
observers which are comoving at a given time. 

Different vacua are related to each other by a Bogoliubov transformation || 

0k{t) = ai k ip k (t) + f3 k p k {t). (19) 

We can from this transformation directly calculate distribution functions and 
expectation values of a field in the |0) vacuum state, seen by an observer corresponding to 
the 0^> vacuum state. Notice that as a consequence of our parametrization equation (JTSJ), 
we can read off the Bogoliubov coefficients for a transformation from the initial vacuum 
state to a comoving state to be a k ,j3 k as defined in equation (JTBJ). 

If we suppose that the particle production is resonant around t , the effective mass 
of tp can be Taylor expanded near to. Furthermore if everything happens approximately 
during one e-fold we will also suppose a ~ a and the equation of motion for the <^-modes 
equation (fT3|) in this neighbourhood is 

(a 3/2 p k ) " + (| + 9 2 \M\t - to) 2 ) (a 3/ V) = 0. (20) 

Now we can rewrite equation (jlHjl in terms of the dimensionless momentum p 2 = 
o~ x k 2 ja\ and the time T 2 = oit — to) 2 ) where a = g\4>o\, 

^ (a 3 /V P ) + {P 2 + T 2 ) (a 3 ^) = 0. (21) 

This is equivalent to a Schrodinger equation for scattering at a (negative) parabolic 
potential and can be solved in terms of parabolic cylinder functions. The incoming 
wave can be approximated as pure positive frequency corresponding to being in the 
vacuum state, while the outgoing state (a k , f3 k ) is found to be (see [T2| for details) 



af = VI + e— P 2 e lCfc , (3° k ut = -ie~* p - ilCfc , (22) 

a = Kgr (l^)4(l-4). (23) 

The number density of ip k modes after particle production is given as ^2] 

a 3 n k (t) = (0 |jV fc | 0) = (0 \a{(t)a k (t)\ o) = |$™T = e"^ 2 = ^1^1. (24) 

If we suppose that the typical time scale of the inflaton field is H^ 1 then O — 0o-^o- 
From the slow-roll result and the limit on the quadrupole we find that if particle 
production has happened at an observationally interesting scale, the exponent above 
becomes 

nH / k/a \ 2 
m { H ) ' 

|| Strictly speaking this is only true for vacua which share the same spatial base |13j . 
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This goes to show that initially the spectrum of 93-modes will be quite blue. Integrating 
equation we find the total number density to be 



n v {t) = j^T-J dk 3 n k = g^-ff^ ( - ) = //„ ( - ) • (25) 



When we try to calculate the two-point correlator (p 2 ) ultraviolet divergences arise. 
These can be dealt with exactly like in normal flat Minkowski-space by subtracting the 
vacuum contribution ^0] and we find 

,0. 1 f d/c/c ^ 

(<?) 



2ir 2 a 3 J uJk 



\(3 k \ 2 + Re (a k Pl exp ^-2i J uj k dt^j 



C nof-) , (26) 



m — gcf) \a Q 
where C = 0.680 is a numerical constant^. 

2-4- The inflaton field 

Energy from the created p particles is draining the kinetic energy of the background 
inflaton field, but in turn the p field reacts back on the fluctuations of the inflaton. The 
inflaton field is a light field, it is frozen on large scales, and it is important to use the 
full equation of motion equation (JHJ) including metric terms for its fluctuations Qk- 

The homogeneous part of the p field is negligible to begin with, because the field 
is overdamped on large scales. If this is the case, and we later disregard oscillations 
between the two fields, i.e. correlators of the form (<pp), then the homogeneous part of 
the (p field is easily seen to remain negligible. A similar analysis of second-order effects 
from a scalar field coupled to the inflaton during preheating |22j found no impact from 
linear perturbations, and indeed showed that the dominant effect is second order in the 
scalar field perturbation. This in turn means that we can disregard any off-diagonal 
terms in the equation of motion for the perturbations of the (f) field. Their equation of 
motion, using equation (jJJJ), is then 



Qk = —3HQ k 



k 2 dV(0) 2Kn 2x 



^ (( SH + i\i + 2^-gX{m-g4>)(< f ?) 



m 2 pl H\\ H 



Qk, (27) 



where we naively have replaced any quadratic terms in <p with the corresponding 
correlators and Qk = 5<t> + (4>/H)ip. 

In a proper treatment one should include other second-order terms in equations ([27)1 . 
(JEJ), © and (fTU]l coming from metric perturbations [TBI EH EHl EH EH], but most probably 
it will only lead to a quantitative change of the produced feature in the power spectrum, 
not a qualitative one. 



C = J p 2 dpc wp (1 — vT-j- c^p 2 sin£fc)/ J p 2 dpe np . The correction C is due to the difference in 
expanding \/3k\ 2 + Re(. . .) correctly (nominator) compared to the approximation used (denominator). 



Observational constraints on particle production during inflation 



8 



2.5. Numerical results 



In order to find the impact on the power spectrum, we carried out a numerical integration 
of the whole system including the two-point correlators such as the one in equation (j2fij) . 
The obvious strategy would be to evaluate the two-point correlators in terms of the 
Bogoliubov coefficients, but this scheme turns out to be numerically unstable. Instead 
we use the idea of [14J and parametrize (p k as 

exp (-i f tOkdt'^ 



[! + /]■ 



28 



Making / the dynamical variable we get the advantage that the two-point correlators 
can be evaluated directly in terms of /. After some algebra, we find all the two point 
correlators: 

k 2 dk M, 



1 



2vr 2 a 3 
1 

4n 2 a 3 
1 

47r 2 a 3 



k 2 dk 

k 2 dk 
2uj k 



2 \f\ 2 + Re(f) 
Re((l + f)f) 



{tA H]w+2Re{f)) 



(29) 
(30) 



\ff + 



(( 



2u k 2 



\f\ 2 + 2Re(f) 



+2Re(f) (u k Im{f) - (1 + Re(f)) ( 

+ lH)lm(f) 



I u k 



+2Im(f)(u k (l + Re(f)) + {^ ! } 
We could instead use tp k directly to compute e.g. (y? 2 ), then 



2vr 2 



k dk 



(31) 



(32) 



2a 3 tv k _ 

It was found numerically that evaluating the correlators directly in terms of (p k is prone 
to numerical instability, because the above integral has to be performed as a discreet 
sum over modes, where even small roundoff or numerical error in ip k for high k leads 
to large errors in the resulting correlator. In effect, using / instead, we only evolve the 
deviation from the vacuum mode. 

Notice that in all of the above expressions the Active contribution from the vacuum 
expectation value has been subtracted, rendering them finite. 

Inserting the definition of / from equation (|2*Hj) into equation ()13|) we find its 
evolution equation 



./ 



2iu k + — J / 



3 /tjfcV 

4 wj 



u k 
2u k 



-H 2 - -H 

4 4 



(33) 



To make any actual integration we have to choose a specific inflationary model by 
specifying the inflaton potential V(4>). For ease of comparison, we have chosen to use 
the same model as Chung et al [TU] : 



V{4>) = \ml4> 2 



m<f> 



pi ■ 



m 



2m 



pi- 



(34) 
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with varying M . To check if the slope of the slow-roll potential had any impact, several 
runs with the exponential potential 

y{<t>) = \ vr KiJ n \i ex P I ol-^— | , a = 0.7034, 
V ' 2 * pl V m pi J (35) 

m = 10~ 6 m p; , m = 2m ph g = 1, 

were carried out. The values of a and are chosen such that and the Hubble 
parameter H are exactly the same as in equation ()34)1 at the moment of particle 
production. This potential is ideally suited to test if the inflaton potential has any 
impact, because it has a substantially higher acceleration of the inflaton. We found that 
the effective step height N e g is virtually unaltered between the two models, confirming 
the resonant nature of the phenomena. 

By starting the simulation at a point so far back in time that all modes of interest 
were way inside the horizon, we used the Bunch-Davies vacuum state equation as 
initial values for (fk and Qk- Furthermore the inflaton was taken to be slow-rolling as 
an initial condition: 



V,^ ((j) start) tr2/n\ _ 



0(0) = fatart, 0(0) = Z vZ \ H (°) = ^V(<Pstart), (36) 

3m; 




/ = 0, / = 0, (38) 

where e(k) is a random variable on the complex unit circle. To check that these values 
had no impact on the further evolution, different starting values <p s tart were chosen for 
different runs, and in general a sta rt was tuned such that at m — g0 o = we have ao = 1. 

To get the proper attractor solution of the background inflaton field, we let the 
simulation run for some time only evolving the homogeneous component of the inflaton 
until the initial solution has relaxed to the correct attractor solution. Afterwards we 
start evolving the (p- and 0-modes. 

Comparison between the analytical approximations and numerical computation 
of the time evolution of the two point correlator (<£> 2 ) has been made with excellent 
agreement (see figure HJ. Notice in figure Q the very steep decline in the correlation 
function just at the beginning. We interpret it as the first short time span, where the 
particles are relativistic. If we look at the approximation equation for the correlation 
function it is singular in the beginning, because we, in making the approximation, have 
neglected the /c-dependence of the energy u^. Obviously we have to take that into 
account when HoAt 1. We then find 

(^)rel ^A^f = ~pF= (-Y 2 - (39) 

2ir 2 a 2 J \/g\<k\ UoJ 

If we now, inspired by the form of u>k, make an ansatz for the correlation function of 
the form 

(f')= ... "° . J -Y\ (40) 
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Figure 1. The number density and the time development of the correlation function, 
which roughly speaking is the same as the backreaction term. Both are compared 
to the analytic approximations, and in general there are very good agreement. The 
boxcar smoothing, was done to emphasise the secular evolution by canceling out the 
oscillations. 



we find, by comparison with equations (}26|) and (}39|) and taking the appropriate limits, 
that 

H At c = f°, ■ (41) 
2yWo| 

Equation (|4f)|) is plotted together with the numerical result in figure H We can interpret 
HoAt c as the number of e-folds where on average the energy of the produced particles is 
stored in relativistic modes. Because of the CMBR limit on the quadrupole, and recalling 
the slow-roll result, in any model of single field inflation we will have HoAt c ~ I0~ 3 g~ 1 / 2 . 

It is important to confirm the findings of jTHj, that the analytic approximation 
for the Bogoliubov coefficients is very precise. The number density is plotted in 
figure ^and shows excellent agreement. Notice that at relatively high fc-values of the 
wavenumber, compared to the horizon size, the number density is still significant. The 
energy density at the moment of particle production per logarithmic interval scales as 
k s nk and therefore the peak energy occurs at values much larger than Hq. This means 
that the process by far and large is local in nature and the impact from the expanding 
background is small. It assures us that even though we have neglected all the metric 
backreaction terms, we still get a result which at least qualitatively agrees with what a 
proper treatment would yield. 

In general our results for the background parameters such as the inflaton velocity 
and the Hubble parameter are very much in line with JU] and we will not repeat their 
analysis here. The reason is that the detailed evolution of the inflaton or the Hubble 
parameter H do not have a direct observable impact. They would only be needed if we 
wanted to do a detailed analytical modelling. 
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2. 6. The curvature perturbation 



The essential quantity to consider is the comoving curvature perturbation 1Z. It is 
found that in the present case there is a qualitative disagreement between a numerical 
evaluation and any analytical prediction, regardless of whether we use the slow-roll |26j, 
the Stewart-Lyth j2Zj or even the semi- numerical method of Leach et al They all 
disagree with our numerical approach because they assume that the impact of any change 
in the fields or fluctuations is local in the sense that it only affects modes which are just 
about to cross the horizon. However, one cannot assume the curvature perturbation 
to be constant on super horizon scales if the entropy perturbation is changing |2T) 129]. 
which is exactly what happens in our case. The homogeneous inflaton field suddenly 
dumps a large amount of energy into ultra relativistic particles. It leads to a suppression 
of the comoving curvature perturbation on all scales larger than the horizon, and the 
formation of a step in the power spectrum (see figure |2J). 



> K l(T a 

p. 



10* I , 

0.001 




k/<toH 



>e 10 
ft 




k/a H„ 



Figure 2. The power spectra from different models. In the right plot the effective 
number of degrees of freedom is so high that a perturbative calculation of the step 
height cannot be trusted, while the left plot shows the monotonic evolution of step 
height for weaker couplings. 



On top of the main feature, the step in the potential, there are high-frequency 
oscillations. As is obvious from the graph the frequency increases with increasing wave 
number k. This suggests a qualitative interpretation, from a mathematical point of 
view, of both the step and the oscillations. The comoving curvature perturbation 1Z is 
related to the inflaton field fluctuation Q k according to equation (J7|) while the power 
spectrum is defined in the usual way 

pi ' 2 =S |K| - <42 » 

To understand features in Vn we can equally well analyse Qk, because the change in H 
is negligible and the change in is only transient (see equation Q). 

While inside the horizon a single perturbation with definite wave number k is 
oscillating with a frequency that is approximately given as k/a (for a massless field). 
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At the moment of particle production we can schematically write the equation for this 
mode Q as 

Q + 3HQ +(K + B(t)\ Q = 0, (43) 

where B(t) is the backreaction term. Now depending on the phase of Q at the exact 
moment of particle production the final amplitude will either be increased or decreased 
[SHj (see second plot in figure EJ). Inside the horizon the phase is highly oscillating 
Q ~ exp(— ik/at) and we get wiggles on the UV part of the spectrum. Modes which 
are far outside the horizon have approximately k/a = and the amplitude changes very 
slowly, therefore they all experience the same kick from B(t) making the step in the 
potential. 

In the short interval, where B l / 2 (t) is greater than the Hubble mass 3/2H, the 
super-horizon modes begin to oscillate again. In the perturbative regime this 'oscillation' 
lasts for such a short time that the modes do not even cross zero a single time, and we 
can compute the impact using simple perturbation theory. On the other hand, if Bit) 
is large enough, zero-crossing occurs and effectively the impact from B(t) is larger than 
the starting amplitude and perturbative methods become meaningless. 

Let us define the effective step height N eS as the quotient between the ultraviolet 

1 12 

and the infrared part of the power spectrum V-^ , which is the main observational 
signature. For small values of the coupling it is possible to make an analytic 
approximation for iV e ff. This is very valuable, because it enables us to connect the 
inner parameters of the theory, i.e. the coupling, the effective degrees of freedom etc, to 
the observable: the effective step height N eS . 

For low ^-values, until the particle production starts, Qk is evolving very slowly 
and can be taken to be constant. Then, suddenly, particle production commences. The 
main term in the equation of motion for Qk equation pTj) is the two point correlation 
function. Define the rescaled mode Q = a 3 ^ 2 Q and use e-folds iV = HAt as time 
variable, then approximately 

^ . » W „, (44) 



dN 2 y H 2 4 / 
and the boundary conditions are, still for low fc-values only, 

Q(0) = Qo, ^f(O) = IQo- (45) 

Unfortunately, inserting the expression equation (pfUj) for the correlation function, we get 
an equation without analytical solutions. But we can exploit the fact that the integrated 
impact of the correlation function is supposed to be small. Then the exact functional 
form is not so important, only the integrated size, and instead of the correlation function 
can we substitute an exponential Ae~ 3N with the correct late time behaviour. The 
equation can now be solved in terms of Bessel functions and plugging in the boundary 
conditions we get the late time behaviour and can evaluate the effective step height 

N cS = 1 - -A + —A 2 + .... (46) 
9 324 v ; 
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A can be determined by the requirement that the integrated size of the perturbation 
has to be the same as for the correlation function: 

oo poo I 

Ae~ 3N dN = g 2 M / dN^J- (47) 
o Jo H 2 



where 7 e ~ 0.577 is Euler's constant and we have evaluated the integral to first order in 
HoAt c . The result is in agreement with the numerics. 

3. Models with a changing EOS 

Clearly the last section shows that the impact on the curvature perturbation on large 
scales is a generic feature. It has to do with the rapid change of the potential energy or, 
alternatively, of the equation of state. We find our result in agreement with analyses of 
similar models, all posessing the property of having a sudden change in the EOS. 

Starobinsky [HIj constructed a phenomenological model, where the potential has a 
break in its slope. Roberts et al considered the so-called false vacuum model with 
quartic potential. This model is characterized by two periods, where first the potential 
energy is dominated by a quartic term in the inflaton, until at a certain moment a 
cosmological constant term takes over the leading role. In the transition period, there 
can be a short suspension of inflation, with related change in the EOS. 

The index of EOS 7 is defined as 

P = ^V 

P <f'/2 + V • 
where we for simplicity have assumed a single-field model. The evolution equation for 
the 0-field is a second-order equation, and it is clear that a change in the EOS can only 
be sourced by a change in the potential energy. In all cases diminishes. Because of 
the equation of motion 

+ 3#0 + = 0, (49) 

this leads to a large second-order derivative and the inflaton enters a stage of fast-roll, 
where approximately 

+ 3#0 « 0. (50) 

After some time the balance between the friction term and the source term is 
reestablished. Looking at equations (0) and (j33J) we can estimate that a significant 
change in power spectrum will only happen when 

^% > 1, (51) 
9H 2 ~ ' K J 

for at least one e-fold. Supposing that V</> during one e-fold is diminished by a factor e, 
the average change during this period is 

^ = -(1 - e)V^ = 3(1 - 6)F0. (52) 
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Using the last two equations we find a bound on e 

e < |. (53) 

It is clear that these are very rough estimates, but it does show that to get a significant 
change in the power spectrum of the fluctuations a rather large change in V$ is required. 
This could either be due to interactions with other fields, as in our case, or due to the 
internal dynamics of the model, as in the two mentioned cases. 



4. Observational data 

4-1- Cosmic microwave background 

The CMB temperature fluctuations are conveniently described in terms of the spherical 
harmonics power spectrum 

Q = (\a lm \ 2 ), (54) 

where 

AT 

—(e,<j>)=Y,aimYi m {OA)- (55) 

lm 

Since Thomson scattering polarizes light, there are additional power spectra coming 
from the polarization anisotropics. The polarization can be divided into a curl-free (E) 
and a curl (B) component, yielding four independent power spectra: Ctj, Ce,i, Cb,i 
and the temperature i?-polarization cross-correlation Cte,i- 

The WMAP experiment has reported data only on Ct,i and Cte,i, as described in 

in u 

We have performed the likelihood analysis using the prescription given by the 
WMAP collaboration which includes the correlation between different Cis pP [Sj . 
Foreground contamination has already been subtracted from their published data. 

In parts of the data analysis we also add other CMB data from the compilation 
by Wang et al 33 which includes data at high I. Altogether this data set has 28 data 
points. 



4-2. Large scale structure 

The 2dF Galaxy Redshift Survey (2dFGRS) [33] has measured the redshifts of more 
than 230 000 galaxies with a median redshift of z m w 0.11. One of the main goals of 
the survey was to measure the galaxy power spectrum on scales up to a few hundred 
Mpc, thus filling in the gap between the small scales covered by earlier galaxy surveys 
and the largest scales where the power spectrum is constrained by observations of the 
CMB. A sample of the size of the 2dFGRS survey allows large-scale structure statistics 
to be measured with very small random errors. An initial estimate of the convolved, 
redshift-space power spectrum of the 2dFGRS has been determined [IT] for a sample of 
160 000 redshifts. On scales 0.02 < k < 0.15/i Mpc -1 the data are robust and the shape 
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of the power spectrum is not affected by redshift-space or nonlinear effects, though the 
amplitude is increased by redshift-space distortions. A potential complication is the 
fact that the galaxy power spectrum may be biased with respect to the matter power 
spectrum, i.e. light does not trace mass exactly at all scales. This is often parametrized 
by introducing a bias factor 

where P s {k) is the power spectrum of the galaxies, and P m (k) is the matter power 
spectrum. Indeed it is well established that on scales less than ~ 10 Mpc different 
galaxy populations exhibit different clustering amplitudes, the so-called morphology- 
density relation (see e.g. [HH]-[38 ). Hierarchical merging scenarios also suggest a more 
complicated picture of biasing as it could be non-linear, scale-dependent and stochastic 
[39]-|44|. However, analysis of the semi-analytic galaxy formation models in [45], as 
well as the simulations in 011 suggest that the biasing is simple and scale-independent 
on large scales k > O.lbh Mpc -1 where the power spectrum is well described by linear 
theory, and we restrict our analysis of the 2dFGRS power spectrum to these scales. Two 
different analyses have demonstrated that the 2dFGRS power spectrum is consistent 
with linear, scale- independent bias jl^l HZj. Thus, the shape of the galaxy power 
spectrum can be used straightforwardly to constrain the shape of the matter power 
spectrum on large scales. 

When looking for steps or other features in the primordial power spectrum using 
the 2dFGRS, one should bear in mind that what is measured is the convolution of the 
true galaxy power spectrum with the 2dFGRS window function [TTj, 

Pconv(fc) oc J P g (k - q)\W k (q)\ 2 d 3 q, (57) 

where W is the window function. In an earlier, simplified attempt to look for steps 
in the primordial power spectrum [IB] it was found that the window function of the 
2dFGRS more or less washed out any features in the primordial power spectrum for 
comoving momenta k < O.lh Mpc -1 . The effect on the class of models considered in 
the present paper is illustrated in figure El for four different values of /cbreak- However, 
combining the 2dFGRS power spectrum with CMB data breaks parameter degeneracies 
that are present if each data set is analysed separately, and therefore a combination of 
large-scale structure and CMB data gives tighter constraints on the primordial power 
spectrum than the CMB alone. 



5. Likelihood analysis 

For calculating the theoretical CMB and matter power spectra we use the publicly 
available CMBFAST package [4*U] . As the set of cosmological parameters we choose Q m , 
the matter density, the curvature parameter, Qf,, the baryon density, H , the Hubble 
parameter, n s , the scalar spectral index of the primordial fluctuation spectrum, r, the 
optical depth to reionization, Q, the normalization of the CMB power spectrum, b, the 
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Figure 3. The effect of convolving matter power spectra with the 2dFGRS window 
function. The black lines are the matter power spectra computed with CMBFAST 
for four different values of fcbreak, the red lines are the corresponding spectra after 
convolution with the window function. The vertical bars are the 2dFGRS power 
spectrum data points. 



Table 1. The different priors on parameters other than N e g and fcbreak used in the 



likelihood analysis. 


Parameter 


Prior 




0.28 ±0.14 (Gaussian) 


h 


0.72 ± 0.08 (Gaussian) 


n b h 2 


0.014-0.040 (top hat) 


n 


0.6-1.4 (top hat) 


T 


0-1 (top hat) 


Q 


free 


b 


free 



bias parameter, and finally the two parameters related to the step height, N e g, and 
location, fcbreak- We restrict the analysis to geometrically flat models Q = Q m + Q\ = 1. 

In principle one might include even more parameters in the analysis, such as r, the 
tensor to scalar ratio of primordial fluctuations. However, r is most likely so close to 
zero that only future high precision experiments may be able to measure it. The same 
is true for other additional parameters. Small deviations from slow-roll during inflation 
can show up as a logarithmic correction to a simple power-law spectrum or 
additional relativistic energy density |HH]-|III| could be present. However, there is no 
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Table 2. Best fit \ 2 and number of degrees of freedom for each of the cases shown in 
figure 0] 



Case 


x 2 


d.o.f. 


X 2 /d-oi. 


(a) 


1466.9 


1388 


1.057 


(b) 


1459.2 


1371 


1.064 


(c) 


1440.9 


1360 


1.059 


(d) 


1432.2 


1343 


1.065 



evidence of any such effect in the present data and therefore we restrict the analysis to 
the 'minimal' standard cosmological model. 

In this full numerical likelihood analysis we use the free parameters discussed above 
with certain priors (see table [TJ determined from cosmological observations other than 
CMB and LSS. In flat models the matter density is restricted by observations of Type 
la supernovae to be Q m = 0.28 ± 0.14 |62j . The current estimated range for Qf,h 2 
from BBN is Q b h 2 = 0.020 ± 0.002 [U3], and finally the HST Hubble key project has 
obtained a constraint on Hq of 72 ± 8 km s -1 Mpc -1 jS]. The actual maginalization 
over parameters other than N e g and fcbreak was performed using a simulated annealing 
procedure 

In figure 0] we show results of the likelihood calculation for various different data 
sets. In panel (a) results are shown for an analysis with all available data, WMAP, 
the Wang et al compilation and the 2dFGRS power spectrum. In panel (b) we show 
results for WMAP + Wang, in panel (c) for WMAP + 2dFGRS, and finally in panel 
(d) for WMAP data only. The figure shows the 68% and 95% exclusion limits for 
the parameters N e s and fcbreak- These two contours correspond to A% 2 = 2.31 and 

= 6.17 respectively Best fit values and number of degrees of freedom for each of 
the four cases are listed in table El 

Several things can be noted from this figure. First, the 2dFGRS power spectrum 
only constrains a step-like feature on scales smaller than k ~ 0.05 — 0.1/i Mpc -1 because 
the window function of the survey washes out any features on scales larger than this. 

Second, the CMB data are mainly sensitive to scales corresponding to I ~ 50-100, 
and only to a lesser extent to features at high /. The main reason for this is that CMB 
data do not directly probe the primordial power spectrum, but rather the primordial 
spectrum convolved with an effective window function. This will be discussed in detail 
in the next subsection. The result of the convolution is that a sharp step in P(k) appears 
as a gradual increase in power over a range in I in the CMB spectrum. It is relatively 
easy to mask this effect by changing other cosmological parameters. However, scales 
at / <> 100 were outside the particle horizon at recombination and therefore much less 
sensitive to changes in most of the cosmological parameters. This is the reason why the 
likelihood analysis shows a very stringent bound in this range, and not in the region 
around the first acoustic peak, as could perhaps have been expected. 

Finally, adding the Wang et al compilation to the WMAP data significantly tightens 
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Figure 4. 68% and 95% confidence exclusion plot of the parameters N e g and fcbreak 
for the four different cases described in the text. The horizontal lines show the range 
in /c-space covered by various data sets. CMB data have been converted from k to I 
using the approximate prescription I ~ 2k/ Hq. 

the bound on iV e fj over most of the range in fc-space, particularly of course at the smaller 
scales where the WMAP data suffer from finite angular resolution. However, in some 
parts of parameter space the constraint actually becomes slightly worse by adding other 
CMB data. This is not inconsistent, it just means that the Wang et al data actually 
favour a small step and therefore push the total likelihood towards higher values of N^. 

In the last panel the effect of the final resolution of the WMAP data becomes 
apparent beyond k ~ 0.05/i Mpc -1 , corresponding to I ~ 450, exactly around the scale 
where the error bars are no longer cosmic variance limited. 

Finally, we note that our analysis is restricted to geometrically flat models. 
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However, since spatial geometry affects the CMB spectrum at scales I <^ 100 it is possible 
that the accuracy with which we have constrained N e g on this scale has been slightly 
overestimated. 



5.1. The CMB window function 

CMB data do not directly measure the underlying primordial power spectrum of 
fluctuations. Rather, they measure the spectrum folded with a transfer function in 
the following sense 

Ci = J ^P(k)AUk), (58) 

where Af(k) is the transfer function taken at the present, r = r , and r is conformal 
time. Following the line-of-sight approach pioneered by Seljak and Zaldariagga j3H| this 
transfer function can be written as 

A(k)= / drS^r^T-To)], (59) 
Jo 

where S is a source function, calculated from the Boltzmann equation, and ji(x) is a 
spherical Bessel function. However, in order to get a very rough idea about the effective 
window function Wi(k) = Af(k)/k of the CMB we approximate S with a constant 
to obtain 

(^rV~0 iorkro^l 
wi(k) oc l\ 1 / (60) 

— for£;T ^/. 

This simple equation shows several things: (a) A feature at some specific 
wavenumber k = k* has the greatest impact on the CMB spectrum at Z* ~ /c*t . For 
a flat, matter dominated universe, To = 2/ Ho, yielding U — 2k*/ Hq- (b) The CMB 
window function is quite broad, and narrow features in P(k) are accordingly difficult to 
detect. 

In the present model the main detectable feature is a step in the power spectrum, 
characterized by the amplitude, N e s, and the location, /vbreak- Adding such a step to an 
otherwise scale-invariant spectrum yields a primordial power spectrum 

P(k)=A + N eS e(k-k hieak ). (61) 

Calculating C\ from equation (|58|). using the window function equation (J60j) then 
gives us 



C, oc < 



f A + N eS 

(I/to) 2 

A | iVeff 
, (I /To) 2 &Lak 



for k hrcak T <; / 

(62) 



+ 71 for ^break^O k, I- 
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Figure 5. Ratio of Ci for model with step to model without step. The full line is the 
precise numerical calculation and the dashed is the approximation in Eq. i|t)3|) . 



Finally, we can define the ratio, R, between the CMB spectrum with the step, and 
one with exactly the same cosmological parameters, but without the step 

iVeff 



R 



1 + 



1 + 



A 



I 



for febreak^O ^ I 

for k hreak T ^ /. 



(63) 



A \ ^brcakTO / 

In figure we show the ratio of two spectra calculated with CMBFAST. The first 
model is a standard fl m = 1 CDM model with h = 0.7 and scale- invariant initial 
spectrum, P(k) = A. The second model is with the same parameters, but with a 
step added, with parameters N eS = A and fcbreak = 0.05 Mpc -1 . The full line shows 
the numerical result, and the dashed line shows the approximation from equation (|63|) . 
Clearly our very simple approximation still captures the essential impact of a step-like 
feature on the CMB spectrum. The reason for choosing a CDM model, rather than the 
observationally preferred ACDM concordance model, is simply that in the former case 
r = 2HQ 1 , making the analytic expression in equation particularly simple. Using 
the ACDM model as the benchmark yields almost exactly the same result. 

The conclusion is that the effective window function wi(k) for CMB smooths out 
the step feature in the underlying primordial power spectrum P(k). This in turn means 
that the ability of CMB data to detect step-like features is significantly degraded. 

A feature which is at first surprising is that the constraint on the step amplitude 
A^efr is strongest at fcbreak ~ 0.01, corresponding to I ~ 90 which is significantly above 
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the scale of the first acoustic peak. The reason is that scales above I ~ 100 were outside 
the particle horizon at recombination and therefore the CMB spectrum at these scales 
is not susceptible to acoustic oscillation effects which can otherwise mask a step in the 
primordial spectrum. At / <> 30 the constraint on N eS becomes gradually weaker simply 
because cosmic variance increases the error bars on C\. Clearly, even though the present 
CMB data are most sensitive to step features at scales corresponding to I ~ 50-100 a 
future experiment like Planck which is cosmic variance limited out to I ;> 2000 would 
be able to put strong constraints on N e g at much higher /, simply because the greater 
measurement accuracy would break the degeneracy between a step and changes in other 
cosmological parameters. 

6. Discussion 

We have presented a calculation of the perturbation to the power spectrum that is 
generated under inflation due to a resonant coupling of the inflaton to a boson. The 
calculation has been made only under the assumption of slow-roll evolution of the 
inflaton before and after the interaction, but no constraints have been put on the specific 
form of the potential. In [TO] the same setup was studied, but with a trilinear coupling 
to a fermion field. Our results for the perturbation to the primordial power spectrum 
are dramatically different to their findings. A step-like feature is formed, and on top 
of this, to the ultraviolet side of the step, there is a highly oscillating transient. For 
weak couplings the size of the step scales like g 5 ^ 2 Af and a suppression of power on large 
scales is observed, while at larger values of the coupling and the dengeneracy this simple 
relation breaks down and suppression or enhancement of the large scale power relative 
to the ultraviolet part can occur (see figure 12}. 

The differences we observe with respect to the results of Chung et al jTUj are easily 
explained from the more limited analysis they performed. Indeed similar numerical 
studies support our conclusion. Easther et al |3U] have investigated the impact from a 
step in the potential of the inflaton and found a very similar feature. They reproduce 
the oscillations seen in our model, but no significant amplification is observed. Leach 
et al undertook a detailed investigation of the false vacuum model [HH] where they 
found a suppression of the infrared part of the spectrum caused by changes in the 
entropy perturbation very much like the one we observe. Starobinsky [31] explored the 
very simplest model, where the inflationary potential has a step, and reached the same 
conclusion. Common for all studies is that the effects seen are due to the self contained 
dynamics of the field, while in our model they are caused by the mutual interaction 
between the inflaton and external fields. This is a major difference, and while the other 
models require fine tuning of the potential, such that the 'feature' occurs in the small 
range accessible to observational cosmology, our model can be naturally supported in 
the context of supersymmetric theories with extra dimensions or superstring theories. 
These theories can have a whole hierarchy of particles with masses comparable to the 
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Planck mass [TOJ E3 EE] • It would then be much more natural to expect at least one 
particle with parameters, such that the produced feature falls into the range open to 
observations JH]- Nonetheless, the presented model should be understood as a realistic 
toy model, while a real model should have some basis in particle physics. 

It would have been desirable to do a full numerical study including higher order met- 
ric terms coming from the gravitational interaction of the fields. We think, though, that 
the overall features of the phenomena are correctly described in this analysis. It can be 
argued, that because the particles are produced relativistically and the peak wave num- 
ber is very large in comparison with the Hubble horizon at that epoch, the process should 
be causal in nature, and therefore not much dependent on the expansion of the universe. 

One feature we did miss in our study, which could be important, is the quantum 
treatment of the inflaton fluctuations and the possibility of cross correlations between 
the fields. The step feature is largely due to the backreaction of the coupled particle 
directly on the inflaton fluctuations, and therefore the cross correlation could be almost 
maximal. This leads to the possibility, that some of the energy contained in the inflaton 
fluctuations oscillates into fluctuations in the other species and is redshifted away [69J. 
This clearly has to be addressed in a future study. 

Using the recently published WMAP data together with other CMB data compiled 
by Wang et al and large scale data from the 2DF collaboration we have found new 
upper limits on a break in the primordial power spectrum. We did the full likelihood 
analysis, and our conclusions are robust: the derived limits changed only slightly while 
varying the amount of data put into the analysis. We find that on scales corresponding 
to k ~ 0. 001-0. 03h Mpc -1 the 2a upper limit on N eS — the relative step height — is ~ 0.3. 
This conclusion is not sensitive to the chosen compilation of data and can essentially 
be derived from the WMAP data only. We also note that the CMB spectrum is mainly 
sensitive to features at scales corresponding to / ~ 50-100 and only to a lesser extent 
features at high /. This is because scales corresponding to I ^ 100 were outside the 
particle horizon at recombination and much less affected by changes in most of the 
cosmological parameters. The constraint is significantly strengthened around / ~ 100 
by adding other data than WMAP. This is indirectly through their tightening of the 
limits on the cosmological parameters translating into less freedom to cancel the effect 
of a break around I ~ 50-100. While observations of the CMB give us a wonderful 
tool to constrain the cosmology of our universe there is no a priori reason to expect a 
smooth almost scale invariant power spectrum. The possibility of inherent noise and 
deviations in the primordial power spectrum gives us a promising prospect to probe 
the earliest history of our universe, but also have the potential of fooling cosmological 
parameter estimation at the per cent level. One should remember this extra uncertainty 
of 'unknown physics', when reading about limits set by new experiments. While so far 
there has been no trace of it, the deviation from full scale invariance and a possible 
bend in the power spectrum, as has been inferred by the WMAP team using their data, 
the 2DF survey and Lyman-alpha forest data, tells us that we are reaching levels of 
precision, where features could begin to crystalize out from the experimental errors. 
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It should be noted that our likelihood analysis for fcbreak and N eS is in fact quite 
general and can be applied to other models which produce a break in the power spectrum. 
Of course the oscillating behaviour of the spectrum close to the break is very specific to 
the present model, but the width of the window functions for both CMB and LSS tends 
to blur the oscillations so that the observational constraints are essentially on the step 
height and position. 

Our result for (fcbreak, A eff ) can therefore be compared for instance to the study of 
|7()j . In that study it was found that the WMAP data show a slight preference for a 
step in the power spectrum at k ~ 2-3 x Mpc -1 . In their analysis it was assumed 
that there was no power above the step scale, which in our language corresponds to 
N e Q — > oo. In fact our WMAP-only analysis shows that the region where N c g — > oo and 
^break — 1.5-2.5 x 10~ 4 Mpc -1 provides a slightly better fit than no break. However, we 
find that this effect is not statistically significant. 

Finally we note that our study only probes the scales measurable by CMB and 
LSS experiments. On much smaller scales limits on primordial black holes provide tight 
constraints on the primordial power spectrum (see for instance [71]) and can therefore 
be expected to provide interesting constraints on step-like features, limiting theories 
predicting the occurence of particle production and phase transitions on a much wider 
range of scales. 
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